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Abstract 

Fermions in D = 4 self-dual Euclidean Taub-NUT space are in- 
vestigated. Dirac-type operators involving Killing- Yano tensors of the 
Taub-NUT geometry are explicitly given showing that they anticom- 
mute with the standard Dirac operator and commute with the Hamil- 
tonian as it is expected. They are connected with the hidden sym- 
metries of the space allowing the construction of a conserved vector 
operator analogous to the Runge-Lenz vector of the Kepler problem. 
This operator is written down pointing out its algebraic properties. 
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1 Introduction 

It is highly non-trivial problem to study the Dirac equation in a curved 
background. One of the most interesting geometries is that of the Euclidean 
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Taub-NUT space since it has not only usual isometries but also it admits 
a hidden symmetry of the Kepler type |], 0. Moreover in the Taub-NUT 
geometry there are four Killing- Yano tensors, connected with three Stackel- 
Killing tensors which represent the components of an analogue to the Runge- 
Lenz vector of the Kepler problem |2], j|, |], |5| . 

In other respects, the Taub-NUT metric is involved in many modern 
studies in physics |], ]2j, || . The Dirac equation in the Kaluza-Klein monopole 
field was studied in the mid eighties |7[ but, to our knowledge, the conserved 
observables of the Dirac theory corresponding to the hidden symmetries of 
the Taub-NUT geometry have not been well discussed in literature. An 
attempt to take into account the Runge-Lenz vector of this geometry was 
done in ||, but their approach is incomplete. 

For this reason the present paper is devoted to the construction of the 
Runge-Lenz vector-operator for the Dirac field in the Taub-NUT background. 
The method we use is to calculate the Dirac-type operators defined with the 
help of the Killing- Yano tensors and then to write down the Runge-Lenz 
operator. 

It is known that space-time isometries give rise to constants of motion 
along geodesies. However not all conserved quantities along geodesies arise 
from isometries. There are prime integrals of motion related to hidden sym- 
metries of the manifold, a manifestation of the existence of the Stackel-Killing 
tensors |], [1(J. A Stackel-Killing tensor of valence r is a tensor k^,,,^ which 
is completely symmetric and which satisfies a generalized Killing equation 

fcOu x .../v;A) = 0. 

There are geometries where the Stackel-Killing tensor can have a certain 
root represented by Killing- Yano tensors. We recall that a tensor / m ... Mr is a 
Killing- Yano tensor of valence r if it is totally antisymmetric and satisfies the 



equation / w ...r^ ; A) = llj. The role of the Killing- Yano tensors becomes 
important in theories with spin as in the cases of the pseudo-classical spinning 
manifolds |I2, |13[ or in quantum theory of fermions [[14] . 



On the other hand, Carter and McLenaghan ]T5[ showed that, in the 
theory of Dirac fermions, for any isometry with Killing vector k^ 1 there is an 
appropriate operator 

x k = -z(^v M - \i»i v K, v ) (i) 

which commutes with the standard Dirac operator defined by 

D s = 7 p V p (2) 
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where V p are the spin covariant derivatives including the spin-connection 
while 7^ are the point-dependent Dirac matrices carrying natural indices. 
The usual Dirac matrices with local (hated) indices will be denoted by 7^. 
We specify that Xk is just the generator of the operator-valued represen- 
tation according which the Dirac field must transform under the isometry 



transformations corresponding to [16] . Moreover, as stated in [15], each 
Killing- Yano tensor of valence 2 produces the non-standard Dirac oper- 
ators of the form 

B/ = -^(/;V,-j7YW (3) 

which anticommutes with the standard Dirac operator D s . We observe that 
the operators ([J) and (||) have specific spin parts. 



2 The Dirac field in Taub-NUT space 

Let us consider the Taub-NUT space and the chart with Cartesian coordi- 
nates x^ (fi, z/, ... = 1, 2, 3, 4) having the line element 

ds 2 = g^dx^dx" = ^dl 2 + V(dx A + Adx 1 ) 2 (4) 

where dl 2 = (dx) 2 = (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 is the Euclidean three-dimen- 
sional line element and A is the gauge field of a monopole. Another chart 
suitable for applications is that of spherical coordinates, x' = (r, 9, 0, x)i 
among them the first three are the spherical coordinates commonly associated 
with the Cartesian space ones, x % (i,j, ... = 1, 2, 3), while x + = ~^x A . The 
real number /1 is the parameter of the theory which enters in the form of 
the function 1/V(r) = 1 + [i/r. The unique non-vanishing component of the 
vector potential in spherical coordinates is = /z(l — cos#). In Cartesian 
coordinates we get divA = and B = rot A = fi^. 

When one uses Cartesian charts in the Taub-NUT geometry, it is useful 
to consider the local frames given by tetrad fields, e(x) and e( defined 
I7|. Their components, have the usual orthonormalization properties and 
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give the components of the metric tensor, g^ v = Va^^u an< ^ 9^ = rf*^ e & el k 
where 77^,3 = 8 & & is the Euclidean metric tensor. Its gauge group, G{rf) = 
SO (4), has the universal covering group G(rf) ~ SU(2) ® SU(2). The four 
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Dirac matrices matrices 7°, that satisfy {7°, 7^} = 2?7 a ^, can be taken as 

0). *-U »)■ < 5 ' 

In addition we consider the matrix 

f = 7W7 4 = ( o 2 _° l2 ) (6) 

which is denoted by 7 in Kaluza-Klein theory explicitly involving the time 
TM. These matrices are self-adjoint and give the covariant basis generators 



of the group G(r)), denoted by S al3 = — i[y a , 7 /3 ]/4. 
The standard Dirac operator is defined as [1T3 

D s = fV, = iW^ ■ P + ~^=^Pa + Z -VWfZ* ■ B (7) 

v V £ 



where V„ are the components of the spin covariant derivatives with local 
indices 

V 4 = tVVp + WVV£ ijk Z*B k , V 4 = -^Pa - % -V\/v£* -B. (8) 



These depend on the momentum operators Pi = —i(di — Aid^) and P4 = 
— i<9 4 , which obey the commutation rules [Pi, Pj] = iSij^B^P^ and [-P,-P 4 ] = 
. The spin matrices giving the spin connection are 

= Si + % -fi , S i = ^e {jk S^ k . (9) 

In our representation of the Dirac matrices, the Hamiltonian operator of 
the massless Dirac field reads [HI 



H = fD s = I W ) . (10) 

Wn / 

This is expressed in terms of the operators n = op — iV~ l P/± and n* = 
up + iV~ 1 P4 where o> = a ■ P involves the Pauli matrices <jj. Denoting by 
the usual covariant derivatives we find that the Klein-Gordon operator 
has the equivalent forms 

A = -V M <TV„ = Vn*ir = VP 2 + ^P 4 2 . (11) 
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3 Conserved observables 



We say that an operator represents a conserved observable if this commutes 
with H. The specific off-diagonal form of H suggested us to introduce the 



Q-operators defined as 113 1 



Q(X) 



H. 




Qc 




(12) 



where Qo = iD s = i^H. These operators allow us to associate Dirac oper- 
ators to the Pauli operators, it, ir*, ap, = a • L and a r = a ■ x/r |T8| . 
The remarkable property of the Q-operators is that if [X, A] = then Q(X) 
commutes with H 1I3J . Other conserved observables can be of diagonal form 
T = diag (T^\ T^) when the condition [T, H] = implies 
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V 



(13) 



and, therefore, [T^\ A] = 0. 

The conserved observables may be found among the operators which com- 
mute or anticommute simultaneously with D s and 7 5 . An important type 
of such operators are the generators of the operator-valued representations 
of the isometry group carried by the spaces of physical fields JO], |TB[. The 
space defined by (]3J) has the isometry group G s = SO (3) ® £4(1), of rota- 
tions of the coordinates x l and translations of x 4 , with the Killing vectors 
ki (i = 1,2,3) and respectively. According to Eq.([l|) the £4(1) generator 
is P4 while other three Killing vectors give the 50(3) generators which are 
the components of the whole angular momentum J = L + S as in the flat 
space-times. The difference is that here the orbital angular momentum is 



x 



x x P — /i-P 4 . 

r 



(14) 



The components of J commute with D s and 7 s and, therefore, they are 
conserved, commuting with H. Moreover, they satisfy the canonical com- 
mutation rules among themselves and with the components of all the other 
vector operators (e.g. coordinates, momenta, spin etc.). 

The first three Killing- Yano tensors of the Taub-NUT space pi 



2e 5 Ae' + e ijk e j A e k 



(15) 
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are rather special since they are covariantly constant. According to Eq.( 
after some algebra, we obtain the Dirac-type operators 



Qi 



(16) 



which anticommute with Q and 7 s , commute with H and obey the N = 4 
super algebra 



{Qa, Qb} = 2S AB H 2 , A, B, ... = 0, 1, 2, 3 

linked to the hyper-Kahler geometry of the Taub-NUT space. 
The fourth Killing- Yano tensor of the Taub-NUT space 



(17) 



./ 



Y 



— f + —rjtijke 3 A e k 



(18) 



is not covariantly constant having the following non- vanishing field strength 
components 

2r 2 

— sin 0. (19) 

Taking into account these covariant derivatives we can calculate its corre- 
sponding Dirac-type operator, Q Y , according to the general rule ([$[), as in- 
dicated in Appendix. The result is 



Q 



Y 









-CTrV- 1 



r 
1— 

/' 



Qo j 



a r 
a r V' x 



(20) 



One can verify that this operator commutes with H and anticommutes with 
Qo and 7 5 . 



4 The Runge-Lenz operator 

The hidden symmetries of the Taub-NUT geometry are encapsulated in the 
non-trivial Stackel-Killing tensors k ifJiU , (i = 1,2,3). They can be expressed 
as symmetrized products of Killing- Yano tensors (|T5| ) and ( |i~8|) Q : 

ki^v — — t(/^ja/ 1 ^ + f\\f lX ^ + -^-{k^k iv + k^yki^). (21) 
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In fact only the product of Killing- Yano tensors /' and / leads to non- 
trivial Stackel-Killing tensors, the last term in the r.h.s. of (^1|) being a 
simple product of Killing vectors. This term is usually added to write the 



Runge-Lenz vector of the scalar (Schrodinger or Klein-Gordon) theory 0, 19 



K = -\v^ v V v = l(PxL-LxP)-^-A + fi-P 4 2 (22) 
2 2 2 r r 

whose components commute with A and obey 

[Li, Kj] = ie ijk K k , [1$, Kj] = «(P 4 2 - &)e ijk L k . (23) 

For the Dirac theory the construction of the Runge-Lenz operators can be 
done using products among the operators Q Y and Qi. This procedure repre- 
sents the only possibility to generate non-trivial conserved operators which 
should be not related to the Hamiltonian, like in Eq.(|l7]), or the Casimir 
operators of the symmetry group of the manifold. 



Using an analogy with the relation (21) from the Taub-NUT geometry. 



let us consider the vector operator iV with the components 

N = ^{Q Y , Q,}- JiP 4 (24) 

which, after simple but tedious calculations based on Eqs. ( |A.8| ) and (|A.9| ) 
from Appendix, can be put in the form 

N = -W{PxJ-JxP-i^xP)W~ 1 --Q{x/r)H 



+ 



nfx(E'xB) P 4 + /i-P 4 2 (25) 



where W = diag (1, \/V) and S* is given by (0). 

Since iV commutes with H, its diagonal blocks, iV^ 1 ) and N^ 2 \ satisfy 
the relations (0). It is straightforward to find that the first block which 
commutes with A, 

A?(D = k + |P 4 , (26) 



contains not only the orbital Runge-Lenz operator ( Eg ) but a spin term too. 
Moreover the components of the operator N satisfy the following commuta- 
tion relations: 

[N u P 4 ] = , [Ni, Jj] = te ijk N k , (27) 
[N u Qo] = , [Ni, Qj] = ie t3k Q k P A (28) 
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[N h Nj] = ie ijk J k F 2 + l -e l3k Q t H , F 2 = P 4 2 - H 2 . (29) 



and 



In order to put the last commutator in a form close to ( pBj ) we can redefine 
the components of the Runge-Lenz operator, /C, as follows 

JQ = N i + ~H~ l (F-P i )Q i (30) 

These operators satisfy the desired commutation rules 

[fCi, H}=0, [/Q, Jj] = ie ijk K. k , (31) 

P 4 ] = , [Ki, Q d ] = is ijk Q k F (32) 

and 

[JCi, fCj] = ie ijk J k F . (33) 

The explicit form of them makes obvious the spin contribution to the 
Runge-Lenz vector, re-confirming the result from pseudo-classical approach 
2(J. We specify that there are no zero modes IJTJJ and, therefore, the 



operator H is invertible such that our definition (j30|) of the Runge-Lenz 
operator is correct. Moreover, as in the scalar case |6|, 0, this can be re- 
scaled in order to recover the familiar dynamical algebras o(4), o(3, 1) or e(3), 
corresponding to different spectral domains of the Kepler-type problems. 

Finally, it is worthy to note that commutation relation (|28|) of the Runge- 
Lenz operator with the standard Dirac operator remains valid even if a mass 
term is included in Eq.(Q). 

Appendix: The operator Q Y 

To evaluate the operator Q Y we start with the first term of (|]) 

where X = a-(xxP) + l = cr L + l + iio r P^ is the operator introduced 
in [21]. Since the Dirac matrices with spherical indices, 7^(2') = e'^(x')7 Q , 
corresponding to the orthogonal coordinates r, 9 and are 

7 r (x') = W^ 1 sin 9 cos (p + sin 9 sin <\> + 7 s cos 9) , (A. 2) 

7 6 *(x / ) = r~ 1 vV(7 1 cos 9 cos <p + 7 2 cos # sin — j 3 sin 9) , (A. 3) 
7 *(ar / ) = (rsin^)" 1 v / l 7 (-7 1 sin0 + 7 2 cos0) (A.4) 
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we find that 

rfrfrfi = -rfrfirf = ...= T~ 2 {sill O^V W f f (A.5) 

is completely antisymmetric in r, 9, (f>. Then, we calculate the second term 
of @, according to ( |I9|) and ( |A.5| ), and embedding the result with (|A.1|) we 

The form ( f2"0"|) and other interesting formulas can be derived using the iden- 
tities [cr r , <jp] = 2ir~ 1 X and 

a P \ = -\a P = % -a ■ (P xL -Lx P) - ^AP 4 • (A.7) 

With their help one can demonstrate that the relations 

2i ( AttX \ 

Q Y Q(X) = -Q(a r X)H + -[ Q _^y AX7f *J_ , (A.8) 



n \ u — v v /\y\- /i — 7= 



OA I -Xtt*X 
Q(X)Q Y = -Q(Xa r )H+-[ Q 1 ] (A.9) 



hold for any 2x2 matrix operator X which commutes with A and V. In 
particular for X = 1 it results [Q Y , H] = 0. 
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